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Semiflexible polymers such as filamentous actin play a vital role in the mechanical behavior of 
cells, yet the basic properties of cross-linked F-actin networks remain poorly understood. To address 
this issue, we have performed numerical studies of the linear response of homogeneous and isotropic 
two-dimensional networks subject to an applied strain at zero temperature. The elastic moduli are 
found to vanish for network densities at a rigidity percolation threshold. For higher densities, two 
regimes are observed: one in which the deformation is predominately affine and the filaments stretch 
and compress; and a second in which bending modes dominate. We identify a dimensionless scalar 
quantity, being a combination of the material length scales, that specifies to which regime a given 
network belongs. A scaling argument is presented that approximately agrees with this crossover 
variable. By a direct geometric measure, we also confirm that the degree of affinity under strain 
correlates with the distinct elastic regimes. We discuss the implications of our findings and suggest 
possible directions for future investigations. 

PACS numbers: 87.16.Ka, 62.20.Dc, 82.35.Pq 



I. INTRODUCTION 

The mechanical stability, response to stress, and lo- 
comotion of eukaryotic cells is largely due to networks 
of biopolymers that collectively form what is known as 
the cytoskeleton 0, 0, 0, 0] • Filamentous actin (F- 
actin), microtubules and other intermediate filamentous 
proteins make up the cytoskeletal network, along with 
a variety of auxiliary proteins that govern such factors 
as cross-linking and filament growth. By understand- 
ing the relation of the individual filament properties and 
dynamically evolving gel microstructure to the rheologi- 
cal/mechanical properties of intracellular structures, one 
will better understand the general framework for cellu- 
lar force generation and transduction 0,10,0. Such 
stress production and sensing underlies such fundamental 
biological processes as cell division, motility 0], and ad- 
hesion 0, 0, 0, 0] . Given the importance of biopoly- 
mer networks in determining the mechanical response of 
cells, there is an obvious interest in understanding the 
properties of such networks at a basic level. Understand- 
ing stress propagation in cells also has implications for 
the interpretation of intracellular microrheology exper- 
iments 0, 0|. However, biopolymers also belong to 
the class of semiflexible polymers, so-called because their 
characteristic bending length (however defined) is com- 
parable to other length scales in the problem, such as the 
contour length or the network mesh size, and thus cannot 
be neglected. 

Such semiflexible polymers pose interesting and fun- 
damental challenges in their own right as polymer ma- 
terials. The understanding of the properties of in- 
dividual semiflexible polymers is quit e hi ghly devel- 
oped [S [H H H M M El IHlll HH; in ad- 
dition, the dynamical and rheological properties of these 



polymers in solution [13, H3, H^, 113, HH H3, H^l have 
largely been elucidated. The remaining problem of deter- 
mining the rheology of permanently cross-linked gels of 
semiflexible polymers, however, has proved quite subtle. 
Related theoretical approaches considered thus far have 
either assumed a simplifying network geometry, such as 
a lattice [34| or a Cayley tree [3{|, or assumed that the 
dominant deformation modes are affine [3r3 | or dominated 
by transverse filament fluctuations and bending (3?i . 

In this paper we study the static mechanical proper- 
ties of random, semiflexible, cross-linked networks in the 
linear response regime with the aim of shedding light 
on the more complex, nonequilibrium cytoskeleton. Our 
approach is deliberately minimalistic: we consider two- 
dimensional, athermal systems with no polydispersity in 
filament properties. Although obviously simplified, this 
restricts the parameter space to a manageable size and al- 
lows for a fuller characterization of the network response. 
The central finding of our work is the existence of quali- 
tatively distinct regimes in the elastic response and local 
deformation in networks, each with characteristic signa- 
tures that should be observable experimentally. 

The basic distinction between stress propagation in 
flexible and semiflexible networks is that in the former 
elastic deformation energy is stored entropically in the 
reduction of the number of chain conformations between 
cross- links [3!j,|4(J while in the latter, the elastic energy is 
stored primarily in the mechanical bending and stretch- 
ing of individual chains. In a flexible, cross-linked mesh 
there is only one microscopic length scale: the mean dis- 
tance between cross- links or entanglements, l c . Since the 
actual identity of a flexible chain is immaterial on scales 
beyond this mean cross-linking distance, chain length 
plays only a very small role 41]. In a semiflexible gel, 
however, chains retain their identity through a cross-link 
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because the tangent vectors on a given chain remain cor- 
related over distances much longer than £, the mesh size. 
Thus, it is reasonable to expect that the elastic proper- 
ties of the network depend on both the mesh size and the 
length of the chains. 

Upon increasing the density of filaments and thus the 
density of cross-links (hereafter we refer to filament den- 
sity only as one quantity determines the other in 2d) the 
system acquires a static shear modulus via a continuous 
phase transition at the rigidity percolation point. This 
corresponds to moving from left to right in the lower 
half of the phase diagram shown in Fig. fl] This critical 
point (solid line in Fig. is at higher density than the 
connectivity percolation point. Since our model ignores 
the entropic elasticity of the network, below this critical 
density associated with rigidity percolation, the material 
has no static shear modulus and may be considered a liq- 
uid. We will discuss the finite temperature implications 
of this zero-temperature phase transition in more detail. 
The elastic properties of the fragile gel (solid) that exists 
just above this critical point are controlled by the physics 
of rigidity percolation. 

By increasing the cross-link density further we en- 
counter a regime over which the elastic deformation of 
the gel is dominated by filament bending and is highly 
nonaffine. This nonaffine regime (NA) is consistent with 
prior predictions by Kroy and Frey [3^ . l38l | . Within the 
NA regime the static shear modulus scales linearly with 
the bending modulus of the individual filaments, k, but 
the elastic moduli are not controlled by properties of the 
rigidity percolation critical point. Most remarkably, how- 
ever, the deformation field under uniformly applied stress 
is highly heterogeneous spatially over long length scales 
comparable even to the system size. We will quantify the 
degree of nonaffinity as a function of length scale and use 
this nonaffinity measure to demonstrate that the degree 
of nonaffinity in the NA regime increases without bound 
as one goes to progressively smaller length scales. Such 
materials are then poorly described by standard contin- 
uum elasticity theory on small length scales. 

By further increasing the filament density, one ap- 
proaches a crossover to a regime of affine deformation 
(A), in which the strain is uniform throughout the sam- 
ple, as the velocity field would be in a simple liquid 
under shear. This cross-over is shown in Fig. ^ by the 
dashed lines above the NA region. In this regime, elastic 
energy is primarily stored in the extension/contraction 
of filaments. Also, in contrast to the NA regime, the 
growth of the degree of nonaffinity saturates as a func- 
tion of decreasing length scale. The elastic response of 
the network is governed primarily by the longitudinal 
compliance of filaments, and the shear modulus can be 
calculated from the combination of this realization and 
the assumption of affine, uniform deformation as shown 
in Refs. [Il HJ EJ M, HI H^. We will show that 
there is one dimensionless parameter A that controls the 
NA — > A crossover. It is set by the ratio of the fila- 
ment length to a combination of parameters describing 
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FIG. 1: A sketch of the expected diagram showing the various 
elastic regimes in terms of molecular weight L and concentra- 
tion c ~ l/l c - The solid line represents the rigidity percola- 
tion transition where rigidity first develops at a macroscopic 
level. This transition is given by L ~ c _1 . The other, dashed 
lines indicate crossovers (not thermodynamic transitions), as 
described in the text. As sketched here, the crossovers be- 
tween nonaffine and affine regimes demonstrate the indepen- 
dent nature of these crossovers from the rigidity percolation 
transition. 



the density of the network and the individual filament 
stiffness. 

In the affine regime, the longitudinal response actu- 
ally can arise from two distinct mechanisms: there are 
two forms of compliance of a semiflexible filament under 
extensional stress, one essentially entropic [3(| and the 
other essentially mechanical ,3Sj . In the first case, the 
compliance relates to the thermally fluctuating filament 
conformation, which, for instance, is straightened-out un- 
der tension. A change in the length of a filament between 
cross-links results not in simple mechanical strain along 
that filament but rather in a reduction of the popula- 
tion of transverse thermal fluctuations along that fila- 
ment thereby reducing the entropy of the filament. This 
reduction results in an elastic restoring force along the 
length of the filament. This is the dominant compliance 
for long enough filament segments (e.g., between cross- 
links). In the second case, the compliance is due to a 
change in the contour length of the filament under ten- 
sion, which, although small, may dominate for short seg- 
ments (e.g., at high concentration). Thus, in general, we 
find two distinct affine (A) regimes, which we refer to as 
entropic (AE) and mechanical (AM). 

Moving still further up and to the left in the diagram 
shown in Fig.^we would eventually reach a regime (not 
shown) in which the filament lengths between cross-links 
are much longer than their thermal persistence length 
and standard rubber elasticity theory would apply. This 
regime is of no experimental importance for the actin 
system. To complete our description of the diagram, we 
note that by increasing the filament concentration in the 
entropic affine regime, one must find a transition from 



3 



entropic to mechanical elasticity (AE — ► AM) within a 
regime of affine deformation. 

Here we confine our attention to two-dimensional per- 
manently cross-linked networks and consider only en- 
thalpic contributions to the elastic moduli of the system. 
In effect we are considering a zero-temperature system, 
except that we account for the extensional modulus of F- 
actin that is principally due to the change in the thermal 
population of transverse thermal fluctuations of a fila- 
ment under extension. We do not expect there to be a 
significant entropic contribution to the free energy com- 
ing from longer length scale filament contour fluctuations 
because of their inherent stiffness. This justifies our ne- 
glect of the sort of entropic contributions to the filament 
free energy that are typically considered in the analy- 
sis of rubber elasticity of flexible polymers. We discuss 
this more below. The low dimensionality of the model 
system, on the other hand, is more significant since the 
essentially straight chains in two dimensions will not in- 
teract sterically under small deformations. In our model, 
only one length scale is required to describe the random 
network, the mean distance between cross-links, l c . In 
three dimensions, however, chains can interact sterically 
and two quantities are needed to fully describe the ran- 
dom network, the density of cross-links and density of 
filaments. 

In our zero-temperature analysis of the system pre- 
sented in this paper, we do not explicitly probe the dif- 
ference between the NA — > AE and the NA — ► AM 
cross-overs. Their difference enters our description of the 
system via a choice made for the form of /it, the exten- 
sion modulus of an individual filament. For physiological 
actin, we expect that the relevant transition will be the 
NA — ► AE. We discuss this further in section Ivl 

The remainder of this paper is arranged as follows. 
In Sec. |n] we define our model system in terms of the 
mechanical properties of individual filaments and the 
manner in which they interconnect to form the network. 
An overview of the simulation method used to find the 
mechanical equilibrium under an imposed strain is also 
given. We describe in Sec. IIIII the rigidity percolation 
transition, at which network rigidity first develops. We 
then describe the crossover from a nonaffine regime above 
the rigidity transition to an affine regime in Sec. IIVI A 
scaling argument is also presented for this crossover. The 
macroscopic mechanical response is also demonstrated to 
be linked to geometric measures of the degree of affine 
deformation at a local level. In Sec. |Vl we show how 
the network response in the affine regime can be either 
essentially thermal or mechanical in nature. In Sec. IVII 
we discuss primarily the experimental implications of our 
results. 



II. THE MODEL 

The bending of semiflexible polymers has been success- 
fully described by the wormlike chain model, in which 



non-zero curvatures induce an energy cost according to 
a bending modulus n. For small curvatures, the Hamil- 
tonian can be written 




where u(s) is the transverse displacement of the filament, 
and s is integrated along the total contour length of the 
filament. Transverse filament dynamics can be inferred 
from this Hamiltonian For finite temperatures, 
can also be used to predict the longitudinal response of an 
isolated filament; however, as T — > the filament buckles, 
preconfiguring a breakdown of the linear response |38| . 
We also consider the response of a filament to compres- 
sion/extensional deformations through the elastic Hamil- 
tonian given by: 

^stretch = ^fljds (^Pj (2) 

where dl /ds gives the relative change in length along the 
filament. (0) is just a Hookean spring response, with 
a stretching modulus /i that is here taken to be inde- 
pendent of k, although they can both be related to the 
cross-sectional radius and elastic properties of individual 
filaments (see Sec. II V (J)) . 

The networks are constructed by the sequential ran- 
dom deposition of monodisperse filaments of length L 
into a two-dimensional shear cell with dimensions W x 
W. Since the position and orientation of filaments are 
uniformly distributed over the allowed ranges, the net- 
works are isotropic and homogeneous when viewed on 
sufficiently large length scales. Each intersection between 
filaments is identified as a cross-link, the mean distance 
between which (as measured along a filament) is denoted 
l c , so that the mean number of cross-links per rod is 
L/l c — 1. Deposition continues until the required cross- 
linking density L/l c has been reached. An example net- 
work geometry is given in Fig. [21 

The system Hamiltonian is found by using discrete ver- 
sions of Q and @ which are linearized with respect 
to filament deflection, ensuring that the macroscopic re- 
sponse is also linear. The detailed procedure is described 
in the Appendix. It is then minimized to find the network 
configuration in mechanical equilibrium. Since entropic 
effects are ignored: we are formally in the T = limit. 
The filaments are coupled at cross-links, which may exert 
arbitrary constraint forces but do not apply constraint 
torques so that the filaments are free to rotate about 
their crossing points. We comment on the validity of this 
assumption in greater detail in Sec. IVII below. Specifi- 
cally, we find that whether or not the physical crosslinks 
are freely-rotating, the mechanical consequences of such 
crosslinks is small for dilute networks, provided that the 
networks are isotropic. Once the displacements of the 
filaments under the applied strain of magnitude 7 have 
been found, the energy per unit area can be calculated, 
which within our linear approximation is equal to 7 2 /2 
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FIG. 2: An example of a network with a cross-link density 
L/l c ~ 29.09 in a shear cell of dimensions W x W and peri- 
odic boundary conditions in both directions. This example is 
small, W = \L\ more typical sizes are W = 5L to 20L. 

times the shear modulus G or the Young's modulus Y, 
for shear and uniaxial strain, respectively 0] ■ Thus the 
elastic moduli can be found for a specific network. The 
procedure is then repeated for different network realiza- 
tions and system sizes until a reliable estimate of the 
modulus is found. See Figs. |31 and 0] for examples of 
solved networks. 

The free parameters in the model are the coefficients 
/i and k and the length scales L and l c . We choose to 
absorb k into a third length scale lb derived from the 
ratio of [i and k, 



Although the simulations assume a constant angular cur- 
vature between nodes, it is possible to assign the phys- 
ical interpretation of the natural length over which a free 
filament bends when differing tangents are imposed at 
each end, as seen by a simple Euler minimization of the 
Hamiltonian. Hence we use the subscript lb, denoting 
'bending length.' Since fi now gives the only energy scale 
in the problem, it scales out, along with one of the length 
scales (say L), and thus we are left with two dimension- 
less control parameters: the filament rigidity lb/L and 
the cross-link density L/l c . Note that there is also a 
fourth length scale, namely the system size W, but all 
of the results presented below are for sufficiently large 
systems that the VF-dependence has vanished. 

We now explore the various deformation regimes of the 
system beginning with the most fragile, sparse shear- 
supporting networks that exist just above the rigidity 
percolation transition. 



FIG. 3: Color online An examples of a low-density network 
with L/l c S3 8.99 in mechanical equilibrium, with filament 
rigidity l\,/L = 0.006. Dangling ends have been removed, 
and the thickness of each line is proportional to the energy 
density, with a minimum thickness so that all rods are visible 
(most lines take this minimum value here). The calibration 
bar shows what proportion of the deformation energy in a 
filament segment is due to stretching or bending. 



III. RIGIDITY PERCOLATION TRANSITION 

For very low cross-link densities L/l c the rods are ei- 
ther isolated or grouped together into small clusters, so 
that there is no connected path between distant parts of 
the system and the elastic moduli vanish. As the den- 
sity of cross-links is increased, there is a conductivity 
percolation transition at L/l c S3 5.42 when a connected 
cluster of infinite size first appears 0] • If there was an 
energy cost for rotation at cross-links, an applied shear 
strain would now induce a stress response and the elas- 
tic moduli would become non-zero |4J|- This is also the 
case when thermal fluctuations generate stresses along 
the filaments |4^. However, for networks with freely- 
rotating cross-links at zero temperature, such as those 
under consideration here, the network is able to deform 
purely by the translation and rotation of filaments. Such 
a floppy mode costs zero energy and thus the elastic 
moduli remain zero. This continues to be the case un- 
til the rigidity percolation transition at a higher density 
L/l c « 5.93 iMlIlEliii, when there are sufficient ex- 
tra constraints that filaments must bend or stretch and 
the moduli become non-zero. 

A full description of the network behavior just above 
the transition has been given elsewhere |50| . so here we 
summarize the results. Just above the rigidity transi- 
tion, both G and Y increase continuously from zero as 
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FIG. 4: Color online The same as Fig. for higher densities 
L/lc « 29.09 (a) and L/l c sa 46.77 (b). For calibration of the 
colors see Fig. [S] 



a power in L/l c , with different prefactors but the same 
exponent /, 
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We have found that / = 3.0±0.2 |50ll51| . consistent with 
the value 3.15 ± 0.2 found independently |44|- It is also 
possible to measure geometric properties of rigid clusters, 
such as their fractal dimension; this has been done using 
the pebble game method, and found to give exponents 
that are similar to that of central force (i.e. Hookean 
spring) percolation on a lattice |46j |. However, such net- 
works cannot support bending, whereas we have found 
that the system Hamiltonian for our model is dominated 



by its bending term near to the transition. We conclude 
that our system is in a different universality class to cen- 
tral force percolation, at least as far as the ideas of uni- 
versality apply to rigidity percolation; indeed, it casts 
doubt on the validity of universality for force percolation 
as a whole, as discussed below. Note that this discrep- 
ancy cannot be due to any form of long range correlation 
in the morphology of the system, since our random net- 
works are constructed in such a way as the ensure geo- 
metric correlations cannot extend beyond the length of 
a single filament. Note also that although cross-links in 
our random networks are connected by filament sections 
of varying lengths, thus producing a broad distribution of 
spring constants which can also destroy universality ac- 
cording to the integral expression in 

Hi 

. our networks do 

not violate this condition. This is simply because there 
is a maximum length L, and hence a minimum spring 
constant (J>/L, between any two cross-links, ensuring a 
low-end cut-off to the distribution of spring constants. 
Similar considerations hold for the bending interaction. 

Of the exponents that are used to characterize the crit- 
ical regime, those measured by the pebble game method 
reflect topological or geometric properties of the grow- 
ing rigid cluster. The exponent / is of a different class 
since it measures the mechanical properties of the fragile 
solid that appears at the critical point. Our observation 
of distinct /'s for two systems (i.e. our simulations and 
central force lattices) that appear to share the same geo- 
metric and topological exponents suggests that, although 
there are large universality classes for the topological ex- 
ponents describing the interactions that produce the ap- 
propriate number of constraints, the scaling of the shear 
modulus admits a larger range of relevant perturbations. 
We suspect that while rigidity itself is a highly nonlocal 
property of the network, the modulus depends critically 
on how stress propagates through particular fragile, low 
density regions of the rigid cluster. Thus the modulus de- 
pends on details of how stress propagates locally through 
perhaps a few cross-links so that the mechanical charac- 
teristics of the filaments and the cross-links become rel- 
evant. 

The possibility of experimentally observing the physics 
associated with the zero-temperature, rigidity percola- 
tion critical point sensitively depend on the size of 
the critical region. As with other strictly zero tempera- 
ture phases transitions, there can be experimental con- 
sequences of the critical point physics only if the system 
can be tuned to pass through the critical region, since the 
critical point itself cannot be explored. At finite temper- 
ature, the network below the rigidity percolation point 
(but above the connectivity, or scalar percolation point) 
has a residual static shear modulus generated by entropic 
tension in the system. One might imagine that one may 
crudely estimate size of the critical regime around the 
rigidity percolation transition by comparing the decaying 
zero-temperature (i.e. mechanical) modulus of the net- 
work above the phase transition G ~ (L/l c — L// C |trans) 
to the entropic modulus below it, G on ~ k^T/l^. Unfor- 
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tunately to make such an estimate one implicitly makes 
assumptions about the, as yet unknown, cross-over ex- 
ponents. Regardless, we speculate that the critical per- 
colation point may indeed have physical implications at 
room temperature. First, due to the significant stiffness 
of the filaments, this residual entropically generated mod- 
ulus should be small and thus the physics of the zero- 
temperature critical point may have experimental rele- 
vance. Secondly, one may interpret the observed differ- 
ence between the numerically extracted and calculated 
scaling exponents z that describes the dependence of A 
on l c /lb as evidence of corrections to mean- field scaling 
due to the proximity of the rigidity percolation critical 
point. This interpretation is strengthened by the fact 
that additional data points at smaller values of l c (i.e. 
higher cross-link density and thus farther from the rigid- 
ity percolation point) but fixed A conform more closely to 
our mean-field scaling exponent. Thus, it is reasonable to 
expect that one may observe phenomena associated with 
rigidity percolation in sparsely cross- linked actin systems. 



IV. ELASTIC REGIMES 

The coarse-grained deformation of a material is nor- 
mally described by the strain field that is defined at all 
spatial points. Both the internal state of stress and the 
density of stored elastic energy (related by a functional 
derivative) are then functions of the symmetrized defor- 
mation tensor. In our model, the underlying microscopic 
description consists entirely of the combination of trans- 
lation, rotation, stretching and bending modes of the fil- 
aments. Of these, only the latter two store elastic energy 
and thereby generate forces in the material. Thus a com- 
plete description of the energetics of a filament encom- 
passing these two modes must lead to a macroscopic, or 
continuum elastic description of the material. 

The state of deformation itself, however, is purely a 
geometric quantity; it can be discussed independently of 
the energetics associated with the deformation of the fil- 
aments themselves. We will characterize the deformation 
field as affine if deformation tensor is spatially uniform 
under uniformly applied strain at the edges of the sample. 
Of course, this strict affine limit is never perfectly real- 
ized within our simulations, but, as shown below, there is 
a broad region of parameter space in which the deforma- 
tion field is approximately affine, in the sense that quan- 
titative measures of the network response asymtote to 
their affine predictions. This entire region shall be called 
"affine". Since in continuum elastic models the stored 
elastic energy depends only on the squares (and possi- 
bly higher even powers) of the spatial gradients of the 
symmetrized deformation tensor, it is clear that uniform 
strain is a global energy minimum of the system consis- 
tent with the uniformly imposed strain at the boundary. 

The spatial homogeneity of the strain field allows one 
to draw a particularly simple connection between the 
elastic properties of the individual elements of the net- 



work and its collective properties. Because under affine 
deformation every filament experiences exactly the same 
deformation, the collective elastic properties of the net- 
work can be calculated by determining energy stored in 
a single filament under the affine deformation and the 
averaging over all orientations of filaments. Note this 
calculation constitutes a mean-field description of elas- 
ticity. 

If the strain field is purely affine (i.e., is uniformly dis- 
tributed throughout the sample) on length scales larger 
than the microscopic length scales (e.g., the distance be- 
tween cross-links), then it can be Taylor expanded to 
give a locally uniform strain in which all elements of the 
strain tensor are constants. It is then straightforward to 
see that the filaments would purely stretch (or compress) 
and the moduli would be independent of filament bend- 
ing coefficient k. Indeed, it is possible to derive exact 
expressions for G and Y in this described below. 

Conversely, nonaffine deformation on microscopic 
lengths arises as a result of filament bending, and hence 
a dependency on n and well as /U. This is what we find 
above the rigidity percolation transition, on increasing 
network concentration or molecular weight. Surprisingly, 
however, this is not restricted to the neighborhood of the 
transition, but constitutes a broad regime of the available 
parameter space. This nonaffine regime is dominated by 
bending, as can be seen by the fact that G and Y are 
independent of /i below. It is important to note that 
continuum elasticity breaks down on length scales over 
which the deformation field is nonaffine. In addition, the 
appearance nonaffine deformations invalidates the sim- 
ple, mean-field calculation of the moduli which assumes 
that every filament undergoes the deformation. Gener- 
ically, the moduli in the nonaffine deformation regime 
will be smaller than their value calculated under the as- 
sumption of affine deformation. Nonaffine deformation 
fields in effect introduce more degrees of freedom since 
the deformation field is nonuniform. Using those extra 
degrees of freedom, the system is able to further lower 
its elastic energy by nonaffine deformations and thereby 
reduce its modulus. Only upon increasing the number of 
mutual constraints in the system, can one constrain the 
system to affine deformations and thereby maximize its 
modulus. 



A. Nonaffine, bending— dominated 

Starting from the most sparse networks just above the 
rigidity percolation transition, we first encounter the non- 
affine regime on increasing L/l c . We find empirically that 
throughout this regime (until the cross-over to affine de- 
formation of one kind or the other - AM or AE) the mod- 
uli of the network are controlled by the bending modes 
of the filaments. This nonaffine response can be distin- 
guished from the scaling regime around the transition, 
which is also nonaffine, by the lack of the diverging length 
scale associated with a continuous phase transition; in 
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the simulations, this corresponds to the independence of 
G and Y on system size W for relatively small W, as 
opposed to the increasingly large W required for conver- 
gence close to the transition. 

The dependence of G on the system parameters in this 
regime can be semi-quantitatively understood as follows. 
Consider what happens when k — > 0. In this limit, fil- 
aments can freely bend and only stretching modes con- 
tribute to the dynamic response. This Hookean case has 
already been investigated by Kellomaki et al. for the 
same random network geometries as considered here |53| . 
They found the striking result that floppy modes exist 
for all densities in the linear response, i.e G = Y = 0. 
In these floppy modes, the filaments will bend at cross- 
links, but without costing energy since k = 0. If k is now 
continuously increased from zero, then there should be a 
range of sufficiently small k in which the angles remain 
unchanged but now incur an energy cost according to Q , 
giving a total energy and hence G that is proportional to 

K. 

To make this idea more specific, suppose that when k 
the angles of filament deflection at cross-links, {59}, 
are distributed with zero mean and variance cr 2 e . By 
assumption, cr| e can only depend on the two geometric 
lengthscales L and l c , or to be more precise on their ratio 
L/l c . From JJJ, the energy at each cross-link is 
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(5) 



The mean number of cross-links per unit area is NL/2l Cl 
where N is the number of filaments per unit area. N can 
be exactly related to L/l c using the expression derived 
in |50|: however, for current purposes it is sufficient to 
use the approximate relation L/l c w (a — 1)/(1 — 2/a), 
where a — 2L 2 N/tt. Thus © summed over the whole 
network gives 



calculate the corresponding modulus G a ffmc . Consider a 
rod of length L lying at an angle 9 to the cc-axis. Under 
a shear 7 xy , this will undergo a relative change in length 
5L/L — 7 xy sin#cos#, and therefore an energy cost 



^stretch = TjA^Txy sin 2 9 cos 2 I 



(7) 



The sin 2 9 cos 2 9 factor reduces to 1 /8 after uniformly av- 
eraging over all 9 G (0, tt). Summing over the network as 
in the nonaffine case gives 
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(8) 



where we have also corrected for dangling ends by renor- 
malising the rod lengths L — > L — 2l c . In the high — 
density limit L/l c — > oo, (JHJ asymtotes to 



^affinc ^ „ , 

16 l c 



(9) 



Since the number of rods per unit area is N ^ l/(Ll c ), 
then the concentration of protein monomers of charac- 
teristic size a is c ~ NL/a ~ (aZ c ) _1 and thus G ~ (ac) a 
with a = 1. For comparison, theories of thermal 3D sys- 
tems predict a = | for calculations based on a tube pic- 
ture [13, , a = Tp for affine scaling relations [3^| , 
and a = 2 for the T = three-dimensional cellular 
foam ; 34]. 

The above calculation can be repeated for an affine uni- 
axial strain 7 yy to give similar expressions for the Young's 
modulus Y, with the sin 2 9 cos 2 9 term in J7J replaced 
with sin 4 9. Since this averages to 3/8, Yaffi ne differs by 
a factor of 3, 



Gu ir^xlM 



(6) 



^affine — 3G a ffi 1] 



(10) 



Therefore plotting GL/fi — (TggLl^/l^ versus lb/L on log- 
log axes will give a straight line of slope 2, as confirmed 
by the simulations below. It is also possible to infer the 
variation of a$g on L jl c from either the simulations or the 
scaling argument presented below, but since this is not an 
easily measurable quantity experimentally, nothing more 
will be said about it here. A bending-dominated re- 
sponse was assumed in the calculations of Frey et al. pisf 
and Joly-Duhamel et al. jH^i although the 2D and 3D 
density dependencies will of course be different. 



B. Affine, stretching— dominated 

Under an affine strain, the network response consists 
purely of stretching modes and it is straightforward to 



Hence the Poisson ratio v — Y/2G—1 for afimely sheared 
networks is v a ffi nc = \ , which should be compared to the 
3D lattice prediction v = | 0] and the 3D Cayley tree 
value v = \ H3- 



C. Scaling argument and crossover between elastic 
regimes 

We now attempt to identify the dominant mode gov- 
erning the deviation from the affine solution. The rele- 
vant length scales for this mode are derived, which, by 
comparing to other lengths in the problem, allow us to es- 
timate when the crossover between stretching-dominated 
and bending-dominated regimes should occur. This pre- 
diction correctly predicts the qualitative trends of the 
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deviation from affinity with the lengths L, l c and lb- 
However, it is not as successful quantitatively as an em- 
pirical scaling law described in Sec. IIV Dl Nonetheless 
we believe it contains the essential physics and therefore 
warrants a full description. 

To proceed, we note that the stretching-only solution 
presented above assumes that the stress is uniform along 
a filament until reaching the dangling end. It is more 
realistic to suppose that it vanishes smoothly. If the rod 
is very long, far from the ends and near the center of the 
rod it is stretched/compressed according to the macro- 
scopic strain 70. We assume that this decreases toward 
zero near the end, over a length lu, so that the reduction 
in stretch/compression energy is of order ^JqI\\ ■ The am- 
plitude of the displacement along this segment, which is 
located near the ends of the rod is of order c? ~ To^|| - This 
deformation, however, clearly comes at the price of de- 
formations of surrounding filaments, which we assume to 
be primarily bending in nature (the dominant constraints 
on this rod will be due to filaments crossing at a large an- 
gle). The typical amplitude of the induced curvature is of 
order d/l\, where l± characterizes the range over which 
the curved region of the crossing filaments extends. This 
represents what can be thought of as a bending correla- 
tion length, and it will be, in general, different from In. 
The latter can also be thought of as a correlation length, 
specifically for the strain variations near free ends. We 
determine these lengths self-consistently, as follows. 

The corresponding total elastic energy contribution 
due to these coupled deformations is of order 



Aid 



-M7o*|| 



h 



(ii) 



where the final ratio of In to l c gives the typical number 
of constraining rods crossing this region of the filament in 
question. In simple physical terms, the rod can reduce its 
total elastic energy by having the strain near the free ends 
deviate from the otherwise affine, imposed strain field. 
In doing so, it results in a bending of other filaments to 
which it is coupled. From this, we expect that the range 
of the typical longitudinal displacement l\\ and transverse 
displacement l± are related by 



/ 3 



/i c - 



(12) 



Of course, the bending of the other filaments will only 
occur because of constraints on them. Otherwise, they 
would simply translate in space. We assume that the 
transverse constraints on these bent filaments to be pri- 
marily due to compression/stretch of the rods which are 
linked to them. These distortions will be governed by 
the same physics as described above. In particular, the 
length scale of the corresponding deformations is of or- 
der Z|| , and they have a typical amplitude of d. Thus, the 
combined curvature and stretch energy is of order 



A£ 2 =/i7 2 2||^ 



7o*|| 
I 2 



(13) 



where, in a similar way to the case above, lj_/l c deter- 
mines the typical number of filaments constraining the 
bent one we focus on here. This determines another re- 
lationship between the optimal bending and stretch cor- 
relation lengths, which can be written 



~ ^b^|| ^c- 



(14) 



Thus, the longitudinal strains of the filaments decay to 
zero over a length of order 



'|| ~ lc 



2/5 



(15) 



while the resulting bending of filaments extends over a 
distance of order 



l± ~ l c 



2/5 



(16) 



The physical implications of Eqs. (|15fl and (|16fl is that 
a length of each filament of order lu experiences non- 
affine deformation and this nonaffine deformation causes 
changes in the local strain field over a zone extending a 
perpendicular distance l± from the ends of that filament. 
Thus when lu becomes comparable to the length of the 
filament, L, the network should deform in a nonaffine 
manner. We will refer to this length along the filament 
contour over which one expects to find nonaffine defor- 
mation as A. 

These results make sense, as increased bending rigid- 
ity can be expected to increase the bending correlation 
length l±, while decreasing the longitudinal correlation 
length lu because of the stiff constraints provided by the 
cross-links. Both lengths, of course, tend to increase 
with decreasing concentration of cross-links, i.e., with 
increasing l c . This scaling analysis assumes, however, 
that l c < In j_. Furthermore, we expect that < l c in 
general. This is because the bending stiffness of a rod 
K ~ Y/r 4 increases with the fourth power of its radius r, 
while /i ~ Y jr 2 increases with the square of the radius, 
where Yf is the Young's modulus of the filament. Thus, 
lh is expected to be of order the rod diameter, which 
must be smaller than the distance between cross-links, 
especially considering the the small volume fractions <p 
of less than 1% in many cases. For rods of radius r in 
three dimensions, we expect that l c ~ a/y/cf>. Thus, lb/l c 
may be in the range of 0.01 to 0.1. This means that we 
expect that In > l±, although both of these lengths are 
of order l c . Thus, the natural dimensionless variable de- 
termining the degree of affinity of the strain is the ratio 
of the filament length to In, the larger of the two lengths 
that characterize the range of nonaffine deformation; i.e. 
when the filaments are very long compared with the effec- 
tively stress-free ends, then most of the rod segments ex- 
perience a stretch/compression deformation determined 
by their orientation and the macroscopic strain. 
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It is possible, however, that l± above may become 
smaller than l c , especially for either very flexible rods 
or for low concentrations. This is unphysical, and we 
expect the bending of constraining filaments above to 
extend only over a length of order l c when lb/l c becomes 
very small. This results in a different scaling of In, given 
by 

h ~ ?JK- (17) 

Although we see no evidence for this scaling, it may be- 
come valid for small enough lb/l c - 
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D. Numerical results for elastic moduli 

We now summarize our numerical results starting first 
at lowest filament densities. Away from the rigidity tran- 
sition, the shear modulus G continues to increase mono- 
tonically with the cross-link density L/l c at a rate that 
only weakly depends on the filament rigidity lb/L, as 
shown in Fig. [5] Indeed, for high densities G approaches 
the affine prediction G a ffi ne which, due to the absence 
of bending modes under affine strain as discussed at the 
beginning of Sec. II VI is independent of lb- The filament 
rigidity lb does, however, influence the crossover to the 
affine solution as will be discussed below. The Young's 
modulus Y and the Poisson ratio v — Y/2G — 1 for a 
range of densities are shown in Fig. |SJ It is apparent that 
v remains close to the affine prediction = i, ex- 

cept possibly near the transition where it takes the value 
0.35 ±0.1 [5(| (here as elsewhere, quoted errors are sin- 
gle standard deviations). Note that, in two dimensions, 
area-preserving deformations have u = 1, so ^ = 0.5 
does not imply incompressibility. For comparison, the 
thermodynamically stable range is — 1 < v < 1 |4^. The 
robustness the affine prediction of the Poisson ratio even 
deep in the NA regime is somewhat surprising and is not 
accounted for in our arguments. 

Varying the ratio lb / L over many orders of magnitude 
at fixed L/l c reveals a new regime in which G cx k, rather 
than G ~ G a fR no oc /i as in the affine regime described 
above. An example is given in Fig. [7J where it can be 
seen that G oc l^ oc k, suggesting that this regime is dom- 
inated by bending modes, a claim that is supported by 
the theoretical considerations presented in Sec. II V Al and 
the work of Frey et al. and Joly-Duhamel et al. |54j . 
We also confirm in the inset to this figure that the regime 
for which G ~ G a ffi no is dominated by stretching modes, 
and that this new regime with G oc k is dominated by 
bending modes, as expected. 

The crossover between the two regimes can be quanti- 
fied by the introduction of a new length A, being a com- 
bination of lb and l c characterized by an exponent z, 



FIG. 5: The normalized shear modulus GL/^i versus the 
cross-link density L/l c for three different bending lengths 
lb/L. The solid line gives the affine solution © and the 
dashed lines adjoining the data points are to guide the eye. 
Here and throughout, errors are no larger than the symbols. 
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FIG. 6: The shear and Young's moduli G and Y for lb/L = 
0.006 against L/l c . The interconnecting lines are to guide the 
eye. (Inset) The Poisson ratio v — Y/2G — 1 against L/l c for 
the same lh/L. 

The ratio L/X can then be used to ascertain which regime 
the network is in, in the sense that A « L corresponds 
to the affine regime, and A > L corresponds to the non- 
affine regime. Note that this is only possible outside the 
neighborhood of the rigidity transition. For densities in 
the approximate range 13 < L/l c < 47, a very good data 
collapse can be found by using Ai = \flfjlb or z = |, 
as demonstrated in Fig. |H1 This empirical relation has 
already been published 51]. However, as clearly evident 
from the figure, it appears to fail for the small number 
of very high density points that we have now been able 
to attain. Conversely, the scaling argument of S ec. IIV 01 
generates the relevant length scale A 2 = l\\ = \Jl 7 c /lb or 
z = |. Although the data does not collapse for this sec- 
ond form, as evident from Fig. it appears to improve 
the overall collapse for larger densities (i.e. further from 
the rigidity transition point). A possible explanation for 
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l b / L L / %2 



FIG. 7: Shear modulus G versus filament rigidity lb/L for 
L/l c ~ 29.09, where G has been scaled to the affine predic- 
tion for this density. The straight line corresponds to the 
bending-dominated regime with Gas, which gives a line of 
slope 2 when plotted on these axes. (Inset) The proportion of 
stretching energy to the total energy for the same networks, 
plotted against the same horizontal axis lb/L. 
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FIG. 8: G/G affine versus L/A with A = f/l*/l b for different 
densities L/l c , showing good collapse except for the highest 
density considered. The enlarged points for L/l c ~ 29.09 
correspond to the same parameters as in Fig. 1121 



this is that A2 is the correct asymptotic length, but does 
not apply for networks of intermediate density, where the 
empirical form Ai is much more successful, perhaps due 
to corrections to scaling from the transition point as dis- 
cussed in Sec. IIIII (This is also the relevant range for 
biological applications). However, our current computa- 
tional resources cannot go to higher densities, and so we 
must leave this question to be resolved at a later date, 
by either improved theory or increased processor speeds. 

Although the goal of this paper is to characterize the 
behavior of semiflexible polymer networks over the whole 
of the parameter space, it is nonetheless instructive to 
also consider parameters corresponding to physiological 
actin networks. The lengths L, l c and l\, for F-actin in 
physiological conditions can be approximated as follows. 



FIG. 9: The same data as Fig. |S| plotted against L/A2 with 
A2 = v/'c/'bi as predicted by the scaling argument in the 
text, showing slight but consistent deviations from collapse 
for this range of L/l c . 



The distance between cross-links has been quoted as l c ~ 
0.1/xm |38j| . which for filament lengths L ss 2pm gives 
a cross-link density L/l c « 20. As already argued in 
Sec. IIV CI we can estimate It, by regarding the filament 
as a solid elastic cylinder with radius r, in which case 
l b ~ r w lOnm. Thus we find that lb/L ~ r/L ~ 10~ 3 , 
which gives L/X » 5. Looking at Fig. (SJ this suggests 
that cytoskeletal networks are in the crossover region. 
Similarly, L/X2 ~ 4 leading to the same basic conclusion. 



E. Spatial correlations 

A further way of probing the degree of affinity of a net- 
work is to consider a suitable spatial correlation function. 
Whatever quantity is chosen, it is clear that there can be 
no fluctuations if the strain is purely affine. This is not 
the case with nonaffine strains, which will induce local- 
ized bending modes that couple to the local geometry of 
the network and thus may fluctuate from one part of the 
network to the next. Thus correlations between fluctua- 
tions should have a longer range when the deformation 
is more nonaffine, qualitatively speaking. The two-point 
correlation function between spatially-varying quantities 
A(x) and B(x) can be generally defined as 

C AB (r) = (A(x)B(x + rn)) - (A(x)) (B(x)) (19) 

where the angled brackets denote averaging over all net- 
work nodes x and direction unit vectors n. Fig.llOlshows 
an example of G ep , where p is the local mass density of 
filaments and e is the energy per unit filament length, re- 
stricted to either stretching or bending energy as shown 
in the key. There is a clear anti-correlation between den- 
sity and both forms of energy at short separations, show- 
ing that the magnitude of deformation is heterogeneously 
distributed throughout the network, being greater in re- 
gions of low mass density and smaller in regions of high 
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mass density. Qualitatively, the network concentrates the 
largest deformations into low mass density regions, thus 
reducing the macroscopic energy cost. This effects both 
bending and stretching modes equally: there is no in- 
creased likelihood of one mode over the other for regions 
of given density, as demonstrated by the collapse of C ep 
for both energy types collapse after normalization, also 
given in this figure. 

The sizes of locally-correlated regions can be inferred 
from the decay of a suitable autocorrelation function, 
such as the combined energy E (stretching plus bending) 
per unit length. Cee{t) is plotted in Fig. ^2f° r differ- 
ent density networks. The trend is for Cee{t) to decay 
more slowly with r for lower L/l c at fixed lb/L, suggest- 
ing larger 'pockets' of non-uniform deformation for lower 
network densities, presumably becoming infinitely large 
at the transition, where the correlation length diverges al- 
gebraically with the known exponent v « 1.17±0.02 |4t| . 
We have been unable to extract a meaningful length scale 
from our Cee{t) data and hence are unable to confirm 
the value of this exponent. 
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FIG. 10: The correlation function C ep (r) between local mass 
density p and energy density e, where e is restricted to either 
stretching or bending energy as shown and p is the length 
of filaments within a radius L/4 of the network point. The 
cross-link density L/l c ~ 21.48, lb/L = 0.006 and both lines 
have been normalized so that \C Ep (r = 0)| = 1. 



F. Measures of affinity 

Intuitively, the degree to which the network deforma- 
tion is or is not affine depends on the length scale on 
which we look. For length scales comparable to the sys- 
tem size, the deformation must appear affine since we 
are imposing an affine strain at the periodic boundaries. 
Only on some smaller length scale might deviations from 
affinity be observed. If the deformation field is nonaffinc 
on length scales corresponding to the microscopic lengths 
L, l c or ^b, then the filaments will 'feel' a locally non- 
uniform strain field and the assumptions leading to the 
prediction of G a ffi no will break down. 

To quantify the degree of affinity at a given length 
scale, consider the infinitesimal change in angle under an 
imposed shear strain between two network nodes sep- 
arated by a distance r. Denote this angle 8, and its 
corresponding affine prediction # a ffinc- Then a suitable 
measure of deviation from affinity on length scales r is 

(A8 2 (r)) = ((8-9 affinc ) 2 ) (20) 

where the angled brackets denote averaging over both 
network points and different network realizations. An 
example of (A9 2 (r)) is given in Fig.^J and clearly shows 
that it monotonically decays with distance, as intuitively 
expected. Also, the deviation from affinity is uniformly 
higher for lower lb/L at the same L/l c , in accord with 
the greater deviation of G from G a ffi n c observed above. 

Although (A8 2 (r)) decreases monotonically with r, the 
decay is slow, almost power-law like over the ranges 
given. This suggests that there is no single 'affinity length 
scale' above which the deformation looks affine, and be- 
low which it does not. However, we can read off the 
degree of affinity at the cross-link length scale r = l c , 



which (after normalising to the strain 7) should by <C 1 
for an affine deformation, and 3> 1 for a nonafnne one. 
This is evident when comparing Fig. ^| to the G/G a ffi nc 
for the same systems in Fig- El :^(A# 2 ('c)) <C 1 does 
indeed correspond to G f=s G a ffi no , and G » Gbcnd for 
^<A0 2 (Z C )>»1. 

The monotonically increasing deviation from affinity 
with decreasing L/X can also be seen using an indepen- 
dent affinity measure, as used by Langer and Liu |55| . 
Consider the displacements {<5x,} of each node i after 
the strain has been applied, relative to their unstrained 
positions. Each of these has a corresponding affine pre- 
diction (5x^ ffine that can be simply computed given the 
node's original position and the type of strain applied. 
Then a scalar measure of the global deviation from affin- 
ity is the root mean square of the difference between the 
measured displacements and their affine values, i.e. 



to 2 = ((8xi - Sxf ffiac ) 2 ) (21) 

where the angled brackets denote averaging over all 
nodes i. In Fig. ^| we plot m/L against L/X, and ob- 
serve the expected monotonic increase of the deviation 
from affinity with decreasing L/X. However, the data 
for different L/l c do not collapse. The problem is that, 
unlike in Figs. |H1 and it is not obvious how m should 
be normalised to give a dimensionlcss quantity; we have 
tried using the lengthscales L, lb, l c and A, but none of 
these generate good collapse. It is likely that some com- 
bination of these lengths will collapse the data, but we 
have been unable to find it empirically, and we have no 
theoretical prediction for this affinity measure. 
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FIG. 11: Autocorrelation of the combined (stretching plus 
bending) energy density E for h/L = 0.006 and the L/l c given 
in the key. The system sizes were W = 15L (L/l c = 13.92), 
10L (18.95 and 21.48), 6±L (31.63) and l\L (77.42). 



FIG. 13: The root mean square deviation of node displace- 
ments from their affine prediction, m/L, plotted against the 
same L/Ai as in Fig. [8] Symbol sizes are larger than errors, 
so the apparent scatter is real. 



1000 



100 



10 



0.1 



0.01 



0.001 



L /X 



= 1.63 
3.50 
7.55 
16.26 
35.03 




0.1 



0.034 = L/L 



FIG. 12: Plot of the affinity measure (A9 2 (r)) normalized to 
the magnitude of the imposed strain 7 against distance r/L, 
for different lb/L. The value of r corresponding to the mean 
distance between cross-links l c is also indicated, as is the solid 
line ^■(A6 |2 (r)) = 1, which separates affine from nonaffine 
networks to with an order of magnitude (the actual crossover 
regime is somewhat broad). In all cases, L/l c ~ 29.1 and the 
system size was W = ¥rL. 



V. THERMAL EFFECTS 

In the above, we have considered a mechanical, purely 
athermal model of networks governed by two microscopic 
energies: (i) the bending of semiflexible filaments, and 
(ii) the longitudinal compliance of these filaments that 
describes their response to compression and stretching 
forces. We have already discussed the role of tempera- 
ture in terms of the formation of a solid via the rigidity 
percolation transition. Now, we examine the transition 
between the AM and AE regimes. 

For a homogeneous filament of Young's modulus Y f, 
the corresponding single-filament parameters k (the 



bending rigidity) and [L (the one-dimensional compres- 
sion/stretch modulus) are determined by Yf and geo- 
metric factors: k ~ Y/a , and [i ~ Y/a 2 , where 2a is the 
filament diameter. From here on, we refer to the later 
(mechanical) modulus as /xm- 

At finite temperature, however, there will be trans- 
verse fluctuations of the filaments that give rise to an 
additional longitudinal compliance. Physically, this com- 
pliance comes from the ability to pull-out the thermal 
fluctuations of the filament, even without any stretch- 
ing of the filament backbone. The corresponding linear 
modulus for a filament segment of length I is [3q 



K£ p 
1? 



a 2 L, 



I 3 



-Mm, 



(22) 



where £ p = n/(kT) is the persistence length. The full 
compliance of such a segment is then l/fi = l/^i M + l/^i T , 
corresponding to an effective linear modulus of 



\ l = 



(23) 



Thus, the thermal compliance dominates for lengths 
larger than %J a 2 £ p , while the segment behaves for all 
practical purposes as a rigid rod with linear modulus fiM 
for lengths smaller than yUiHp. 

Thus, there appear to be two distinct regimes within 
the affine regime: for higher concentrations, specifi- 
cally for l c < yf a 2 £ p , the longitudinal compliance of 
the filaments is governed by the mechanical compres- 
sion/stretching of filament segments and the modulus is 
given by Eq. J3J); this is the AM regime. At lower concen- 
trations, specifically for l c > yf a 2 £ p , the single-filament 
compliance is dominated by thermal fluctuations and 



G 



ttk£ p 
T6F' 



(24) 



This is the AE regime. For the networks under discus- 
sion, which are described by a single variable l c that both 
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represents the spacing of filaments and distance between 
constraints or cross-links, the boundary between these 
two affine regimes is simply determined by concentra- 
tion, which is naturally measured as l// c . For actin, we 
estimate the characteristic length \faFl v to be of order 
100 nm. Thus, only when the distance between cross- 
links is less than a distance of order 100 nm will the bulk 
response of the network depend on the purely mechanical 
extension of actin filaments. 

Experimentally, these two affine regimes should be dis- 
tinguished by their scaling dependencies of the linear 
shear moduli G on various parameters. A clearer, quali- 
tative distinction, however, should be seen in their non- 
linear behavior. Specifically, in the thermal regime, the 
maximum strain (either at which the network yields or 
first exhibits non-linear behavior) is expected to decrease 
with increasing concentration of polymer or cross-links 
[3^| . This is because of the limited extent of thermal com- 
pliance, which decreases for shorter filament segments 
that appear more straight. This is in contrast with me- 
chanical networks where the nonlinearities are governed 
by geometry (e.g., connectivity and orientation of fila- 
ments) . This would suggest a concentration- independent 
maximum or characteristic strain, as is seen in some col- 
loidal gels [561 ]. 

Moreover, the actual form of the force-extension rela- 
tion for a semiflexible polymer in the limit of segment 
lengths I <C lp takes on a universal form, depending 
only on a characteristic extension I 2 jl v and characteris- 
tic force tt 2 k/1 2 . This force-extension relation predicts a 
universal strain stiffening of semiflexible gels in the affine 
regime, in contrast with a nonlinear modulus in the me- 
chanical affine regime that will be much more dependent 
in the type of filament [57j . 



VI. IMPLICATIONS AND DISCUSSION 

Starting with solutions having no (zero-frequency) 
elastic behavior, as the filament concentration c/, molec- 
ular weight (filament length) L, or density of cross-links 
is increased, there is a point where macroscopic elastic 
behavior is first observed. This is the rigidity percolation 
transition, and it occurs for a fixed value of L/l c , where 
l/l c is a particularly convenient measure of filament con- 
centration, as it represents the line density (length per 
volume) in our 2-d system, apart from a factor of order 
unity. 

We have shown that there are three distinct elastic be- 
haviors of semiflexible networks above the rigidity tran- 
sition: (i) when either the molecular weight (filament 
length) or concentration is low, a nonaffine regime is 
expected in which the modulus is determined at a mi- 
croscopic level by filament bending (transverse compli- 
ance) [34], |38| ; (ii) as either the molecular weight or con- 
centration increases, a crossover is expected to an elas- 
tic regime in which the deformations are affine (uniform 
strain) and in which the elastic response at the large scale 



is governed by the thermal/entropic longitudinal compli- 
ance of filament segments; (iii) at still higher concentra- 
tions or cross-link densities, this single-filament compli- 
ance becomes dominated by the mechanical compliance 
of bare filament stretching and compression. 

The crossover between (i) and (ii) is given by a fixed 
value of L/X of order 10, where A is a microscopic length 
characterizing the range of nonaffine deformation along 
a filament backbone. We expect this length to be of 
order l c , the distance between cross-links. But, it should 
also depend on the filament stiffness through the length 
lb = \J k/li. In fact, it can only depend on the two 
lengths lb and l c in our networks. Thus, we expect that 
A = l c (l c /lb) z - We have presented a scaling argument 
that shows this for z = 2/5, while we find empirically that 
z 21 1/3 for biologically relevant densities. The boundary 
between nonaffine and affine regimes is thus given by L ~ 
A. In the mechanically-dominated regime, lb = y/ k/ lim, 
while in the thermal regime, 

ib = v^~J-Jr-- ( 25 > 

y a z i v LiM 

Thus, in the mechanical regime, the boundary is given 
by L r~j ll +z , while in the thermal regime, the boundary 
is given by L ~ l\~ z ^ 2 . In either case, we see that this 
crossover has a different functional dependence on con- 
centration, demonstrating once again that the physics 
of this crossover is distinct from the rigidity percolation 
transition. We show a sketch of the expected diagram of 
the various regimes depending on L and c/ in Fig.^ The 
boundary between mechanically-dominated and thermal 
regimes is simply given by l c ~ y/a 2 £ p , as we have noted 
above. 

We can make several additional observations concern- 
ing the behavior of real networks, based on our simple 
model. First, we note the strong dependence of the shear 
modulus on the cross-link density, as illustrated in Eg. 1241 
and already noted for 3D affine networks in Ref . |36J . In 
fact, in the 2D networks presented here, we have made no 
distinction between the mesh size (or typical separation 
of neighboring filaments) £ and the cross-link separation 
l c . We observe from this strong dependence of the modu- 
lus on l c , which is independent of filament concentration 
in 3D, that the modulus of semi- flexible gels can be var- 
ied significantly (in fact, by orders of magnitude) with 
changes only in cross-link density at the same filament 
concentration. This is very different from the situation 
for flexible polymer gels, and may well be important for 
cells, in that the mechanical properties can be tuned by 
local variations in the densities or binding constants of 
various actin binding proteins. 

Furthermore, as we have shown the modulus 

becomes a very strong function of concentration (which, 
again, will translate to cross-link density on 3D net- 
works) in the nonaffine regime. In the nonaffine regime, 
the modulus can vary by several orders of magnitude with 
respect to the stiffer affine gels. As this nonaffine regime 
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is expected for just a few (specifically, of order 10 or 
fewer) cross-links along a single filament, it may also be 
possible that the cell can reduce its stiffness significantly, 
and even fluidize, by decreasing the number of cross- 
links per filament or the filament length. In addition, by 
using its proximity to the NA — ► AE cross-over, the 
cell can tune its nonlinear mechanical properties. In the 
AE regime the cytoskeletal network should be strongly 
strain-stiffening due to the nonlinear extensional prop- 
erties of individual filaments In the NA regime, 
there should be a much larger linear regime since the 
bending modes of the filaments, which dominate the de- 
formations in the NA regime have a much larger linear 
response regime. Finally, we speculate that in the affinc 
regime, the mechanical properties of the cell should be 
insensitive to the details of the cytoskeletal microstruc- 
ture; in the AE regime the mean-field character of the 
network enforced by the large ratio of L to A suggests 
that local effects of cross-linker type or network topol- 
ogy to self-average. On the other hand, within the NA 
regime, the cellular mechanical properties may be quite 
sensitive to such local network modifications. 

In this model, we have assumed freely-rotating 
crosslinks. In the case of actin networks, however, it is 
well know that many associated proteins can bind actin 
filaments at either preferred or fixed angles This 
can have two distinct effects: one geometric, and the 
other mechanical. One the one hand, the model we have 
described is only for isotropic networks. Thus, if actin 
crosslinking results in an anisotropic network {e.g., with 
oriented bundles), then one cannot describe such a sys- 
tem with the model presented here. 

If, on the other hand, the networks remain isotropic, 
but with rigid bond angles between filaments, then we 
expect to see additional rigidity of the networks as a re- 
sult. The size of this effect can be estimated for the 
affine regime (AE and AM, in fact). We find that the 
relative contribution to the network elastic modulus due 
to such crosslinks is only of order ib/^o which is expected 
to be small for realistic networks of actin, both in vitro 
and in vivo. In simple terms, this is simply due to the 
very large lever arm that a filament segment of length 
l c (say, of order 1 micrometer) between crosslinks has 
over a small a small actin binding protein of size a few 
nanometers. More precisely, this can be seen by noting 
that in shearing two filaments that cross at a finite angle 
in the shear plane, a fixed bond angle between the fila- 
ments will give rise to a distortion (i.e., nonaffine) of the 
resulting filament conformations within a region near 
the crosslink. (This length corresponds to the range over 
which finite bending occurs, e.g., when a finite bending 
moment is imposed at a filament end.) The angle of this 
bend will be at most of order the macroscopic strain e in 
the affine regime. This results in a bending elastic energy 
of order ne 2 /l^ per crosslink, compared with the longitu- 
dinal elastic energy of order /ie 2 Z c per segment between 
crosslinks. Noting that I 2 = n/fi, we find that the latter 
term (which corresponds to the freely-rotating crosslink 



case) is larger than the former by a factor of order lb/l c - 
For simple mechanical networks at T — 0, as we con- 
sider in most of this paper, is of order the molecular 
diameter, which is much smaller than the distance be- 
tween overlapping filaments, let alone crosslinks in any 
real actin network. In the in vitro networks that have 
been studied, we expect this ratio to be no larger than 
at most a few percent. In the case of networks at finite 
temperature, as we discuss in Sec. E] lh ~ yWl^vi f° r 
which the ratio above is of order \Jl c /l p . By definition, 
this is smaller than 1 for semiflexible networks. Thus, in 
any case, the corrections to the affine elastic moduli due 
to possible fixed-angle crosslinks (in isotropic networks) 
is expected to be small. 

In the related studies by Wilhelm and Frey ^4}, who 
also consider the T — mechanical properties of net- 
works such as ours, the authors looked at both fixed- 
angle and freely-rotating crosslinks. They found that the 
rigidity percolation transitions occurred at somewhat dif- 
ferent values of concentration for fixed-angle and freely- 
rotating bonds. But, they report that very similar be- 
havior was observed for the two cases above the critical 
points. Specifically, they found no statistically significant 
difference in the dependence of the shear moduli with 
concentration in the non-affine elastic regime. Thus, it 
would appear that no substantial differences due to the 
mechanics of crosslinks can be expected in either affine 
nor nonaffine regimes, at least for the relatively sparse 
isotropic networks that actin forms. 
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APPENDIX 

The filaments are deposited into the shear cell as al- 
ready described in Sec. [H] This is internally represented 
by the set of points {x.;} consisting of all cross-links 
and midpoints between cross-links (the midpoints are in- 
cluded so that the first bending mode between any two 
cross-links is represented). Relative motion between the 
Xj contributes to the system Hamiltonian according to 
discrete versions of Q and (J2J. A change in separation 
from Iq to Iq + SI between any two adjacently connected 
points incurs an energy cost 

(SI \ ^ 
Yj '° (26) 
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In addition to this, a non-zero angle 59 between the vec- 
tors Xj — Xj_i and x i+1 — Xj , where X,_i, x$ and x i+1 are 
consecutive adjacent points on the same filament, con- 
tributes 

5H bead = 1 (27) 

where V is the mean of the lengths to either side of the 
central point, i.e. V — ^(|xj — Xj_i| + |x i+1 — Xj|). cross- 
linked filaments are coupled by imposing the same X; 
at intersections, but there is no energy cost for relative 
angles between filaments: cross-links can freely rotate. 
Each contribution (|26(1 . (|27J) is linearised with respect 



to changes in the x^ and summed to create the system 
Hamiltonian 7i({xi}). Either a uniaxial or shear strain 
7 is applied to the system throu gh t he periodic bound- 
aries in a Lees-Edwards manner |58|. The Hamiltonian 
7i({xi}) is then minimized with respect to the {xj} by 
the conjugate gradient method [H^ . Two optimizations 
are included. The Hessian matrix Aij = d 2 Tl/ dxidxj is 
preconditioned by M , where M has the same diagonal 
2x2 matrices of A but is zero elsewhere. Furthermore, 
cross-links that lie within a given small distance, typ- 
ically w 10 _3 L, are coalesced. This improves the con- 
ditioning of A and hence the speed of convergence con- 
siderably, while producing only minimal change in the 
measured quantities, except precisely at the transition. 
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